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(Dated: February 7, 2008) 

We present measurements of branching fractions and CP- violating asymmetries in decays of B 
mesons to two-body final states containing a K'^ . The resuhs are based on a data sample of 
approximately 88 million T(4S') —> BB decays collected with the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAG. We measure B{B+ K°tt+) = (22.3 ± 1.7 ± 1.1) x_10"'', 
B{B° K°TT°) = (11.4 ± 1.7 ± 0.8) X 10-^ B{B+ ^ K^K+) < 2.5 x 10-^ and S(B" -> K°K°) < 
1.8 X 10~®, where the first uncertainty is statistical and the second is systematic, and the upper 
limits are at the 90% confidence level. In addition, the following CP-violating asymmetries have been 
measured; AcpiB+ ^ K"tt+) = -0.05 ± 0.08 ± 0.01 and Acp{B° ^ K"-k") = 0.03 ± 0.36 ± 0.11. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



The decays of B mesons into charmless hadronic fi- 
nal states provide important information for the study 
of CP violation. In particular, the study of the two- 
body decays B ^ tttt, B ^ Ktt, and B KK pro- 
vides crucial ingredients for measuring or constraining 
the values of the angles a and 7, defined by the ratios 
of various elements of the Cabibbo-Kobayashi-Maskawa 
quark-mixing matrix 1]: a = aig[-V^^V*f^/V^^V*^] and 
7 = arg [— VudKTb/KrfKl]- t'^is paper, we present mea- 
surements of the branching fractions for B meson decays 
to the charmless two-body final states A'°7r+, K^K^^ 
K^TT^, and K^K^ (unless explicitly stated otherwise, 
charge conjugate decay modes are assumed throughout 
this paper and branching fractions are averaged accord- 
ingly). For the B+ K°n+ and B° K°tt° modes we 
also report measurements of the direct CP asymmetries 
in the decay rates: 



Acp = —7= 



r{B^ Pj^TiB^ f) 



T{B^f)+riB^f) 



(1) 



Measurement of the rates and charge asymmetries for 
B Ktt decays can be used to establish direct CP 
violation and to constrain the angle 7 The decay 
B~^ —> K'^tt'^ is dominated by the b —> s penguin process 
and in the Standard Model (SM) is expected to have Acp 
close to zero (< 1%) Thus, observation of a sizable 
charge asymmetry could be an indication of non-SM con- 
tributions to the penguin loop 0, ^ • The B KK de- 
cays are characterized by penguin and VF-exchange pro- 
cesses similar to those in B^ — > tt+tt" and can be used 
01 to determine the angle a from the measurement of the 
time-dependent asymmetries in B^ — > tt+tt". Measure- 
ments of the branching fractions for these decay modes 
also provide important information regarding rescat- 
tering processes. 

The measurements presented in this paper are based on 
data collected with the BABAR detector at the PEP-II 
asymmetric-energy e^e~ collider located at the Stan- 
ford Linear Accelerator Center. The sample consists of 
87.9 ± 1.0 million BB pairs produced at the T{AS) res- 
onance ("on-resonance"), which corresponds to an inte- 
grated luminosity of about 81 fb^^ . An additional 9 fb~^ 
of data recorded at an e+e~ center-of-mass (CM) en- 
ergy approximately 40MeV below the T{AS) resonance 



( "off- resonance" ) are used for background studies. 

The BABAR detector is described in detail in Ref. 
Charged-particle (track) momenta are measured in a 
tracking system consisting of a five-layer, double-sided 
silicon vertex detector and a 40-layer drift chamber 
(DCH), which operate in a solenoidal magnetic field of 
l.brmT. Particles are identified as pions or kaons based 
on the Cherenkov angle measured with a detector of in- 
ternally reflected Cherenkov light (DIRC). The direction 
and energy of photons are determined from the energy de- 
posits in a segmented CsI(Tl) electromagnetic calorime- 
ter (EMC). 

Hadronic events are selected on the basis of charged- 
particle multiplicity and event topology. We reconstruct 
B-meson candidates decaying to K'^X, where X refers to 
7r+, tt", K^ or K^ . The K^ and 7r° candidates are recon- 
structed in the modes K^ K'^ tt+tt" and tt" 77, 
respectively. The following selection criteria are applied 
to the candidate B-decay products. 

Charged tracks are required to be within the tracking 
fiducial volume and to have at least 12 DCH hits and 
a minimum transverse momentum of 0.1 GeV/c. Tracks 
that are not K'^ decay products are also required to orig- 
inate from the interaction point, to be associated with at 
least six Cherenkov photons in the DIRC and to have a 
Cherenkov angle within 4ct of the expected value for a 
pion or kaon. 

Candidate K'^ mesons are reconstructed from pairs of 
oppositely charged tracks that form a vertex with tt+tt" 
invariant mass within 3.5 cr of the nominal K'^ mass and 
measured proper decaytime greater than five times its 
uncertainty. 

Candidate tt^ mesons are formed from pairs of pho- 
tons having invariant mass within 3 a of the nominal tt^ 
mass, where the resolution is about 8 MeV/c^ for the can- 
didates of interest. Photon candidates are required to 
not be matched to a track, to have an energy of at least 
30MeV, and to have the lateral shower shape expected 
for a photon. The 7r° candidates are then kinematically 
fit with their mass constrained to the nominal 7r° mass. 

The B-meson candidate is characterized by two nearly 
independent kinematic variables, the energy-substituted 

mass mEs = \/ is/2 + pi ■ psf / Ef - P% ^^'^ energy 
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difference AE' = Eg — \/s/2, where the subscripts i and 
B refer to the initial e+e~ system and the B candidate, 
respectively, the asterisk denotes the T(45') rest frame, 
and ^/s is the total CM energy. The pion mass is as- 
signed to all charged particles in calculating E^. For 
B'^ — > K^K'^ and i?° K'^tt^ candidates, we require 
|A£;| < 0.11 GeV and |A£;| < 0.15 GeV, respectively. 
For 5+ — > K'^h^ candidates, where h refers to tt or K , 
we require -0.115 < A£: < 0.075 GeV. The interval is 
asymmetric in order to select both B^ K'^tt^ and 
K^K^ decays with nearly 100% efficiency. The 
Ai? distribution is peaked near zero for the modes with 
no charged kaons and shifted on average —45 MeV for 
5+ K^K^ decays due to the pion mass being used 
for the charged B daughter in the calculation. The dis- 
tribution of TTiEs peaks near the B mass for all modes, 
and we require 5.20 < tties < 5.29 GeV/c^. 

Simulated events 0, off-resonance data, and events 
in on-resonance toes and Ai? sideband regions are used 
to study backgrounds. The contribution from other B- 
meson decays is found to be negligible. The primary 
background is from random combinations of tracks and 
neutral clusters produced in the e^e~ — > qq events, where 
q = u, d, s, or c. In the CM frame, this background is 
characterized by its jet structure, in contrast to the more 
uniformly distributed decays of the B mesons produced 
in the T(4S') decays. We exploit this topological dif- 
ference to suppress such background. We require that 
the angle 9g between the sphericity axes of the B candi- 
date and of the remaining particles in the event, in the 
CM frame, satisfies |cos0J| < 0.8. We also construct a 
Fisher discriminant J- given by an optimized linear com- 
bination of J2iPi a'^d J2iP* cos^ 9* lOj, where p* is the 
momentum of particle i and 9* is the angle between its 
momentum and the _B-candidate thrust axis, both cal- 
culated in the CM frame. The shapes of JF for signal 
and background events are included as probability den- 
sity functions (PDFs) in the fits described below. 

Signal yields and charge asymmetries are determined 
from unbinncd extended maximum likelihood fits. The 
extended likelihood for a sample of N K^X candidates 



N 



£ = exp - ^ 



^N,V^{xJ■d,) 



(2) 



where Vi{xj]di) is the probability for a signal or back- 
ground category j, given by a product of PDFs for the 
measured variables Xj of candidate j. The parame- 
ters di determine the expected distributions of measured 
variables in each category and n.i are the yields deter- 
mined from the fit. We perform separate fits for each 
of the three samples of B candidates: B^ — > K^ir^^ 
B" K°K°, and B+ K^h+ {h+ = tt+ or K+). 
For the two neutral B samples there are two categories, 
signal and background, and the yield in each category is 



obtained by maximizing the likelihood. For these fits the 
probability coefficients Ni are the yields (i.e., Ni = rii). 
The charged B decays, B^ K'-^h^, are fit simulta- 
neously with two signal categories, B'^ K'^tt'^ and 
5+ K^K^, and two corresponding background cat- 
egories. In addition, the probability coefficient for each 
category i is given by iV^ — ni{l — qjAi), where nt is the 
total yield, summed over charge states, Ai is the charge 
asymmetry, and qj is the measured charge of the given B 
candidate. The total yields and charge asymmetries are 
determined by maximizing C. 

The independent input variables to the fit xj for a 
given event j are toesj Ai?, and JF. For the fit to 
the B+ K°h+ sample we include the normalized 
Cherenkov residuals {9c^ 9^) /(Jq^ and (^9c — 9^) /crg^, 
where 9c is the measured Cherenkov angle of the primary 
daughter ag^ is its error and 911 (9^) is the expected 
Cherenkov angle for a pion (kaon). The quantities ag^, 
9^, and 9^^ are measured separately for negatively and 
positively charged pions and kaons from a control sample 
of _D" — > K^TT^ originating from decays. 

The parameterizations of the PDFs are determined 
from a combination of data and simulated events. The 
signal TOES PDFs for B+ K°h+ and B^ K°K° 
are derived from fully reconstructed B^ — > D^n^ decays 
and are Gaussian. For B^ K^n^ , simulated signal 
events are employed and the toes PDF is modeled as a 
Gaussian distribution with a low-side power-law tail. We 
use an empirical threshold function |ll| to describe the 
background TOes PDFs. The single shape parameter of 
this function is a free parameter in the B^ K^h^ and 
B^ — > if'^Tr" fits, where the event sample is sufficiently 
large. For the B^ K^K'^ fit this shape parameter is 
determined from on-resonance events in Ai? sidebands. 

The T distribution for signal is modeled as a Gaus- 
sian function with an asymmetric width |l2l |. where the 
parameters are determined from simulated events. For 
background, it is modeled as a sum of two Gaussian 
functions with parameters determined from on-resonance 
events in toes sidebands. 

The signal Ai? PDFs are derived from simulated events 
and are parameterized as a sum of two Gaussian func- 
tions for the modes B+ K°h+ and B° K°K°, and 
as a Gaussian distribution with a low-side power-law tail 
for B'^ K^Ti'^ . The Ai? distribution for background is 
modeled as a second-order polynomial whose parameters 
are determined from on-resonance events in toes side- 
bands. The normalized Cherenkov-angle residuals are 
modeled as a sum of two Gaussian functions. 

The results of the maximum likelihood fits are sum- 
marized in Table HJ The K'^K^ final state is an equal 
admixture of K'^K^ and K^Kl. We therefore assume a 
50% probability for the K°K° to decay as K°K° in com- 
puting the B'^ K^K^ branching fraction. We also use 
the current world averages 13] for B{K2; tt+tt") and 
yB(7r° 77) in computing the branching fractions given 
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TABLE I: Summary of results for numbers of selected K X candidates A^, total detection efficiencies e, fitted signal yields Ns, 
statistical significances S, charge- averaged branching fractions B, charge asymmetries Acp, and 90% confidence-level (C.L.) 
allowed asymmetry intervals. The efficiencies include the branching fractions for intermediate states {K" ■K'^n~ and 

TT 77). Branching fractions are calculated assuming equal rates for r(4S') ^ B°B° and B+B". Upper limits for the K°K+ 
and K^K'^ branching fractions correspond to the 90% C.L. and the central values are given in parentheses. 



Mode 


N 


e{%) 


Ns 


S{a) 


S(10'^) 


Acp 


Acp (90% C.L.) 


K%+ 
K°K+ 


8047 


13.0 ±0.3 
12.8 ±0.3 


255 ± 20j:" 


22 
1.7 


22.3 ± 1.7 ± 1.1 
< 2.5(l.l±0.75i5;;i8) 


-0.05 ±0.08 ±0.01 


[-0.18,0.08] 




2668 


8.6 ±0.5 


86 ± 13 ± 3 


12 


11.4 ± 1.7 ±0.8 


0.03 ±0.36 ±0.11 


[-0.59,0.65] 




754 


8.7 ±0.3 


4.3l^;i ± 1-1 


1.0 


< 1.8(0.61°;^ ±0.1) 







in Table n 

FigureQlshows distributions of mEs and AE for _B+ 
K°Tr~^ and -B" ^ K^tt° candidates after selecting on 
probability ratios to enhance the signal purity. The solid 
curves represent the fit projections after having corrected 
for the efficiency of the additional selection. The efficien- 
cies for these ttt-es (A-E) selection criteria are 70% (93%) 
and 65% (98%) for the K°tt+ and K°tt° states, respec- 
tively, as determined from simulated signal events. 

Signal significance is defined as the square root of the 
difference between — 2 ln£ for the best fit and for the null- 
signal hypothesis. The upper limit on the signal yield for 
a given mode i is defined as the value of n"' for which 

Jq' -Cmaxdrii/ Lmaxdui = 0.9, whcrc £max is the like- 
lihood as a function of rii , maximized with respect to the 
remaining fit parameters. Branching fraction upper lim- 
its are then calculated by increasing the signal yield up- 
per limit and reducing the efficiency by their respective 
systematic uncertainties. 

For the i^T^Tr" mode, which is a CP eigenstate, 

we measure the time-integrated CP asymmetry by deter- 
mining whether the other B meson in the event decayed 
as a i?° or B'^ (flavor tag). The tagging algorithm is de- 
scribed in Ref. |Q]. The measured asymmetry .Amcas is 
given by Acp/{1 + xf), where Xd = 0.755 ± 0.015 13] is 
the B^ mixing parameter. The dilution of the CP asym- 
metry by the factor 1/(1 + z^) is due to the effect of 
B^-B^ mixing in the time evolution of the coherent B^ 
B^ system. 

Systematic uncertainties in the signal yields arise pri- 
marily from imperfect knowledge of the PDF shapes. 
Such systematic errors are evaluated either by varying 
the PDF parameters by their measured (1 cr) uncertain- 
ties or by substituting alternative PDFs from indepen- 
dent control samples. The dominant systematic uncer- 
tainty of this type is that associated with the signal 
Fisher discriminant for both B+ K'^tt^ (±7.1 events) 
and B^ K^tt^ (±1.4 events). Also contributing to the 
systematic uncertainties in the branching fraction mea- 
surements are the uncertainties in the Kl and tt*^ effi- 




5.225 



5.275 



0.1 
AE (GeV) 



FIG. 1: Distributions of thes and for (a,b) B+ K°-k+ 
and (c,d) B'^ candidates that satisfy an optimized re- 

quirement on the signal probability, based on all the variables 
except the one being plotted. The solid curves are projec- 
tions of the fit, while the dashed curves show the background 
contribution. 



ciencies, which are about 3% and 5%, respectively. The 
systematic uncertainties in the charge asymmetries are 
evaluated by adding in quadrature the contributions from 
PDF variations and the upper limit on intrinsic charge 
bias in the detector (±0.01). For the measurement of 
Acp in the decay i?° K^n^ , there is an additional 
contribution of ±0.07 due to uncertainties in the tagging 
efficiencies and mistag fractions. 

In summary, we have measured the branching fractions 
and CP-violating charge asymmetries for 5+ K^tt^ 
and B^ K'^tt'^. No evidence of direct CP violation 
has been observed. We have also searched for the decays 
K°K+ and B° K°K° and set upper limits on 
their branching fractions at 2.5 x 10~^ and 1.8 x 10~^, 
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respectively, at the 90% C.L. The branching fraction 
measurements reported here are consistent with previous 
measurements of the same quantities 0, 0, 0|, but 
have nearly twice the statistical precision. Our measured 
K'^tt'^ charge asymmetry is of the same statistical 
precision and consistent with the value recently reported 
[l8j by the Belle collaboration. All of the aforementioned 
results supersede our previous measurements apart 
from the -B" K'-^tt'^ charge asymmetry, which has not 
previously been measured. 
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We present measurements of branching frac- 
tions and CP-violating asymmetries in decays of 
B mesons to two-body final states containing a 
K^. The results are based on a data sample of 
approximately 88 million T{AS) BB decays 
collected with the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAC. We mea- 
sure B{B+ K^n+) = (22.3 ± 1.7 ± 1.1) x 10"^ 
B{B° K%^) = (11.4±1.7±0.8)xl0-^ 8{B+ 
K^K+) < 2.5 X 10^^ and B{B^ K^K^) < 1.8 x 
10^^, where the first uncertainty is statistical and 
the second is systematic, and the upper limits are 
at the 90% confidence level. In addition, the fol- 
lowing CP-violating asymmetries have been mea- 
sured: AcpiB+ K%+) = -0.05 ± 0.08 ± 0.01 
and AcpiB° ^ K%^) = 0.03 ± 0.36 ± 0.11. 
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